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1. Introduction 



In recent years, research on the AdS/CFT correspondence [1] and its apphcations has expanded 
significantly, leading to substantial new results. At present, holographic descriptions of non- 
perturbative phenomena range from applications to condensed matter systems (e.g superconductiv- 
ity, superfluidity, quantum Hall effect) to applications in high energy physics and to the quark-gluon 
plasma (e.g. confinement/deconfinement phase transitions, chiral symmetry breaking, elliptic flow, 
hadronization). Despite the remarkable insights gained by studying holographic gauge theories, 
the application of the correspondence to phenomenologically relevant gauge theories continues to 
remain a challenge. A way to further improve the phenomenological relevance of the AdS/CFT 
correspondence is to explore further phenomena of universal nature. The results of such studies are 
expected to capture aspects of the quahtative behaviour of realistic gauge theories (such as QCD) , 
which do not possess holographically dual supergravity backgrounds. 

An important example in this class of phenomena is the effect of mass generation and spontaneous 
chiral symmetry breaking in the presence of an external magnetic held. This effect has been shown 
to be model-independent and therefore insensitive to the microscopic physics underlying the low 
energy effective theory. Its essence is the dimensional reduction D ^ D — 2, {3+1 — )■ 1+1) in 
the dynamics of fermion pairing in a magnetic field. The enhanced infrared divergences in lower 
dimensions suggests that the dynamics of fermion pairing is governed by the lowest Landau level, 
which hints at the universal nature of the phenomenon. Furthermore, one can show that the 
{D — 2) -dimensional dynamics of the Landau level favours the condensation of the fermion pairs. 
This effect is known as Magnetic Catalysis. Magnetic Catalysis has been demonstrated in various 
(1+2)- and (l+3)-dimensional field theories [2]-[7] using conventional field theoretical methods. 
The holographic approach to this effect has been initiated in [8], where the (l+3)-dimensional 
holographic gauge theory dual to the D3/D7-brane intersection has been considered^. Additional 
holographic studies of magnetic catalysis at finite temperature or chemical potential for both (1+3)- 
and (l+2)-dimensional systems have been performed in [9]-[23]. 

In the holographic description of magnetic catalysis, the flavour degrees of freedom arc introduced 
by an additional stack of flavour D-branes. The most understood and widely applied regime of 
this approach is when the flavour branes are in the probe limit and their backreaction to the 
ambient supersymmetric background is neglected [24]. On the held theory side, this corresponds 
to the "quenched" approximation in which the number of flavour degrees of freedom, Nf, is much 
smaller than the number of colour degrees of freedom, Nc. In terms of Feynman diagrams this 
implies ignoring the contribution from internal quark loops (windows) in the planar diagrams of the 
corresponding large N^. expansion. At present all holographic studies of magnetic catalysis are in 
the quenched approximation. An obvious question is to ask how corrections due to internal quark 
loops would influence the effect of magnetic catalysis. Our goal is to provide such an estimate in a 
perturbative expansion in the ratio between the number of ffavour and colour degrees of freedom, 

^For a comprehensive review we refer the reader to [10]. 
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For obtaining an unquenched holographic description of magnetic catalysis, we have to take into 
account the backreaction of the flavour branes to the supergravity background sourced by the colour 
branes. Ideally such a background would describe localized branes, however for technical reasons 
this is a difficult task even in the supersymmetric case. To circumvent these difficulties, the flavour 
D-brancs may be distributed along the compact directions of the supergravity background. This 
procedure is called smearing^. The smearing restores a significant part of the global symmetry of 
the geometry and hence simplifies the corresponding Einstein equations. 

Supersymmetric backgrounds of smeared massless flavour D-brancs have been constructed in [27], 
and for flavours with flnite bare mass in [28]. In the case of massive flavours, these backgrounds 
display a hollow cavity in the bulk of the geometry, where the supergravity solution is sourced 
solely by the colour branes. The radius of this cavity is related to the bare mass of the fundamental 
flavors. In the limit of vanishing bare mass the cavity shrinks to the radius of the compact part 
of the geometry and the supergravity background has an essential singularity at the origin of the 
non-compact part of the geometry. In both cases the dilaton held diverges at large radial distances. 
This corresponds to the Landau pole that the dual field theory develops in the UV, due to its 
positive beta function /3 oc Nf/Nf.. 

We can imagine that a non-supersymmetric background interpolating between two supersym- 
metric backgrounds, corresponding to massless flavors in the UV and massive flavors in the IR, 
would describe dynamical mass generation. The radius of the hollow cavity then corresponds to 
the dynamically generated constituent mass of the fundamental flavors. A promising framework for 
the construction of such a geometry was developed in [29], where the ten-dimensional black-hole 
solution dual to the non-conformal plasma of flavoured J\f = 4 supersymmetric Yang-Mills theory is 
presented^. The authors outline the smearing procedure, derive the corresponding equations of mo- 
tion and present a perturbative solution for general massless non-supersymmetric flavour D7-brane 
embeddings. 

The flrst steps towards unquenching the holographic description of magnetic catalysis have been 
undertaken in [33], where the authors, following the approach of [29], construct a perturbative 
non-sup crsymmetric background with a non- vanishing B-ficId, which corresponds to an external 
magnetic fleld coupled to the fundamental degrees of freedom of the dual gauge theory. In the case 
of massless fundamental fields and sufficiently strong magnetic field, the supergravity background is 
unstable, suggesting that the theory undergoes a phase transition to a stable phase with dynamically 
generated mass for the matter fields. 

In section 2 of this paper we complete the studies initiated in [33] constructing a perturbative 
non-supersymmetric background with a non-trivial i?-field, for massive fiavour fields. Our solution 
has a hollow cavity in the bulk of the geometry where it is very similar to the supergravity dual of a 
non-commutative supersymmetric Yang- Mills theory, [34] & [35] . The difference is in the presence 

^For a detailed review on the smearing see [25] while for other solutions employing this technique see [26]. 
^AU the hydrodynamic transport coefficients of the model were analyzed in [30], while the addition of a finite 
baryon density was presented in [31]. For a review on unquenching the Quark Gluon Plasma see [32]. 
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of a squashed S^, instead of a non-squashed one, which breaks the supersymmetry. As suggested 
above, the radius of this cavity, r^, is related to the dynamically generated mass of the fundamental 
fields. 

For radial distances greater than the solution is characterized by a non-vanishing density 
for the smeared D7-brane charge. At sufficiently large radial distance our solution approaches 
the supesymmetric one, constructed in [28]. Following the prescription of [29], we introduce an 
additional large radial parameter ^ (corresponding to a finite UV cutoff ), at which we match 
our solution to the supersymmetric one. Furthermore, we identify the value of the i?-field at as 
the magnetic field of the dual gauge theory, if* = B{r^). 

For radial distances greater than r^, the supergravity background is well approximated by the 
non-perturbative supersymmetric background [28] . This enables us to relate non-perturbatively the 
UV parameters of the theory, namely the finite cutoff h.uv r* and the energy scale corresponding 
to the landau pole of the theory Alp oc tlp, where tlf is the radial distance at which the dilaton 
field diverges. 

Our supergravity construction has the following renormalization group fiow interpretation: 

At the energy scale set by the finite cutoff {Auv c< r») the dual gauge theory is a commutative J\f = 
4 Supersymmetric Yang-Mills theory coupled to Nf fiavours oi Af — 2 hyper mult iplet fundamental 
fields. The fundamental hypermultiplets are coupled to a constant external magnetic field H^, 
which breaks the supersymmetry. Decreasing the energy scale the Yang-Mills theory becomes non- 
commutative and the parameter of non-commutativity (roughly the non-trivial part of the B-Qeld) 
is proportional to the ratio Nf/Nc. 

At energy scales of the order of the physical mass of the fundamental fields, Mg (roughly Mq ~ r^) 
the flavour fields decouple (the D7-brane charge density vanishes). At lower energy scales (inside the 
cavity, r < r^) the dual gauge theory is a pure (only adjoint degrees of freedom) non-commutative 
Yang-Mills theory. To leading order the parameter of non-commutativity in the plane perpendicular 
to the magnetic field scales as ©^^ °^ 

Note that in the bare Lagrangian of the dual gauge theory the external magnetic field couples 
explicitly only to the fundamental degrees of freedom. Therefore the non-commutativity of the 
adjoint degrees of freedom cannot be captured by the quenched approximation and is one of the 
novel results of our analysis. 

Finally, in section 3 of this work we apply our construction to study the effect of magnetic 
catalysis. We develop an appropriate renormalization scheme and compute the free energy and the 
fundamental condensate of the holographically dual gauge theory as a function of the bare mass 
of the fundamental degrees of freedom. Our studies show that to leading order in a perturbative 
expansion in the ratio Nf/Nc, the free energy and fundamental condensate of the theory agree with 
the results obtained in the quenched approximation. Furthermore, at next order in Nf/Nc we show 
that the effect of magnetic catalysis is enhanced and the contribution to the condensate of the 
theory from internal fundamental loops runs logarithmically with the finite cutoff A.uv- 
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2. Constructing the Background 



In the present section we will construct the supergravity background necessary for the holographic 
study of the phenomenon of magnetic catalysis. The field theory duals are realized on the inter- 
section between a set of N(, colour D3-branes and a set of Nf, homogeneously smeared, flavour 
D7-branes, with an additional coupling between the fundamental fields and an external magnetic 
field. The colour D3-branes are placed at the tip of a Calabi-Yau (CY) cone over a Sasaki-Einstein 
manifold X5, where the latter can be expressed as a U{1) fiber bundle over a four-dimensional 
Kahler-Einstein base (KE). The flavour D7-branes extend along the radial direction, wrap a sub- 
manifold X3 of X5 and smear homogeneously over the transverse space [36, 37]. 

2.1 Ansatz &; smeeiring of the flavours 

To take into account the contribution from internal fundamental loops in the Veneziano limit of 
the dual gauge theory, we need to consider the backreaction of the fiavour branes. Ideally the 
corresponding supergravity solution would describe localized D7-branes which break the global 
symmetry of the internal subspace of the supergravity background from S0{6) down to 5*0(4) x 
SO{2) (when is S^). However, even in the supersymmetric case (no external magnetic field), 
this is a very difficult task. One way to circumnavigate the technical difficulties is to construct a 
solution with smeared D7-branes. 

In general the smearing procedure involves distributing the branes at different locations in the 
transverse space subspace and consider a course grained approximation in which the sum over all 
individual embeddings becomes an integral over a distribution of branes. It is somewhat analogous 
to the smearing of point-like charges in electrostatics in 1-1-3 dimensions to obtain: one, two or three 
dimensional charge densities. 

In our case the smearing is performed in such a way that the isometrics of the fibered Kahler- 
Einstein space are kept unbroken, allowing for an ansatz where all the unknown functions just 
depend on a single radial coordinate. Based on this assumption we adopt the following ansatz for 
the metric 

dslf^ = h-'^ [-dt^ + dxl + b{dxl + dxD] + /i^ [b^S^F^da^ + S^dslp2 + F\dT + Acp^)^] , (2.1) 
where the CP^ metric is given by 

dslp2 = \dX^+\ COS^ ^ (d^2 ^ g- ^^^2) ^ 1 ^Qg2 X g- j^2 ^ ^ g^^y ^ 

4 4 2 4 2 2 

1 y 

Acp2 = - cos^ ^(o?V + cos 9dcp) . (2.2) 

The range of the angles is < (%, ^) < tt, < r < 27r, < ^ < An. The ansatz for the NSNS 
and the RR field strengths is given by 

B2 = Hdx"^ Adx^ , C2 = JdtA dx^ , 
F5^Q,{1 + *)e{S'), Fi^Qfp{a){dT + Acp2), F^ ^ dC2 + B2 A F, , (2.3) 
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where £{85) is the volume element of the internal space'' and QcQf arc related to the number of 
different colours and flavours in the following way 

QcVol^Q „ _ 4QfVoliX,) 

^ Vol{X,)g. ■ ^^-'^ 

In our case X5 = and the volume of the three sphere is 27r^. All the functions that appear in 
the ansatz, h,b, S, F,^, J ^H, depend on the radial variable a only. In the convention we follow, 
S &iF have dimensions of length, p,b,h, J ^ H are dimensionless and a has a dimension of length"^. 
Furthermore, a = 00 at the origin and decreases to a* at the boundary. 

The function b in the ansatz for the metric reflects the breaking of the 5*0(1, 3) Lorentz symmetry 
down to 5*0(1, 1) X 5*0(2). The function p{a) in and F(3), determines the distribution of the 
brane embeddings and has a characteristic asymptotic behavior. In fact p{a) encodes the bare mass, 
niq and the fundamental condensate of the dual gauge theory. It vanishes at energy scales smaller 
than the quarks' mass while it asymptotes to 1 in the UV. This leads to the formation of a spherical 
cavity inside the bulk of geometry. The radius of this cavity sets the energy scale related to the 
physical mass of the quark, Mq^ . 

We can understand better the structure of the distribution function p{a) if we consider a repre- 
sentative D7-brane embedding, the so called "flducial embedding". The flducial embedding is an 
auxiliary D7-brane embedding, which probes the backreacted geometry, its shape determines the 
distribution function p{(j). 

To leading order in a multipole expansion we can describe the proflle of the flducial embedding 
by functions, which depend only on the radial (holographic) coordinate a. The flducial embedding 
relevant for our study wraps an internal three-cycle parameterized by 9,(p,ip, extends along a and 
sits at a flxed value of r. To obtain the distribution function p{a) we smear the flducial embedding 
by acting with the symmetries of the internal space (see appendix of [29]). In this way we obtain 

p{a) = cos^ ^ . (2.5) 

To clarify the relationship between the distribution function and the proflle of the flducial embedding 
(2.5), in flgure 1 we have presented plots of a family of D7-brane embeddings obtained by acting on 
the flducial one with a discrete subgroup of the symmetry group of the internal space. Figures 1(a) 
and 1(b) represent probe D7-brane embeddings terminating at the same minimal radial distance 
Tq. The flrst one is a supersymmetric embedding (vanishing magnetic fleld), while the second one 
is a non-supcrsymmetric one (flnite magnetic fleld) and exhibits mass generation. Figures 1(c) and 
1(d) represent the smearing of the corresponding probe D7-brane embeddings. One can clearly see 
the formation of a spherical cavity in the bulk of the geometry. On can also compare the radial 
distributions corresponding to supersymmetric and non-supersymmetric flducial embeddings^. 



^with / e{S^) = Vol(5'5) = 7r3 

^We remind the reader that beyond the quenched approximation the dual field theory has a positive ^-function. 
As a result the mass of the fundamental fields runs with the energy scale and even in the supersymmetric case the 
physical mass differs from the bare mass m^. 

^Notc that we have taken advantage of the fact that to leading order in the ratio Nf /Nc the profile of the fiducial 
embeddings is well approximated by the profile of the probe embeddings. 
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(a) A couple of supersymmetric flat embeddings of (b) A couple of non-supersymmetric mass generat- 

mass TTig — Mq ^ Vq. Note that in the quenched ing embeddings of vanishing bare mass niq = and 

approximation there is no mass generation. physical mass Mq ^ Vq. 




(c) Smeared supersymmetric embeddings. (d) Smeared mass-generating embeddings. 



Figure 1: A visualization of the smearing procedure. 

In the next subsection following the strategy initiated in [29] we will describe the whole system in 
terms of a one-dimensional effective action, from which all the equations of motion can be produced. 
In principle the same system of equations can be explicitly derived also from a ten-dimensional point 
of view after inserting the above ansatz in the ten-dimensional equations of motion plus the Bianchi 
identities. 
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The equations of motion and the corresponding Bianchi identities for the NSNS and RR fields 
are relatively easy to obtain. However, writing down the ten-dimensional Einstein equations is a 
difficult task. The crucial step is to obtain an effective ten-dimensional expression for the smeared 
the DBI action, in the case of massive non-supersymmetric probes. While in the massive super- 
symmetric case such a construction is possible using calibrated geometry (see e.g. [25, 38]) in the 
non-supersymmetric case this is a non-trivial task. To circumvent this difficulty we derive the equa- 
tions of motion from a one dimensional effective action. The key point is that even if we were able 
to obtain a simple ten dimensional term for the smeared DBI action we would still adopt the ansatz 
(2.1)-(2.3) for the supergravity fields, which is equivalent to reducing the action to one dimension. 
Therefore wc could first reduce the eight dimensional DBI action of every individual fiavour branc 
to one dimension and then sum them. In this context the statement that the flavour branes are 
smeared is equivalent to the statement that their reduced one dimensional actions are equivalent. 
In other words, the smearing of the DBI action is equivalent to reducing the eight dimensional DBI 
action of a fiducial fiavour brane embedding to one dimension and multiplying the result by the 
number of flavour branes Nf. 

2.2 Effective action &c the equations of motion 

The action for the Type IIB supergravity plus the contribution from the Nf D7-branes in the 
Einstein frame is 

S = Sub + Sfi , (2.6) 
where the relevant terms of the Sub action are 



■>IIB 



1 1 

23! 



(3) 



-V / C4A//3AF3, (2.7) 



and the action for the flavour D7-branes takes the usual DBI-I-WZ form 



S 



fi 



Nf 



(i^a;e*V-det(G + e-*/2j^) 



Cs + Cq A B2 



(2.8) 



with T = B + 27ra'F. In those expressions B denotes a non-constant magnetic field, F the world- 
volume gauge field and the hat refers to the puUback of the quantities, along the worldvolume 
directions of the D7-brane. The gravitational constant and D7-brane tension, in terms of string 
parameters, are 

1 Ty 1 



2^20 



The first step in our analysis is deriving the equation of motion for the fiducial embedding, which 
follows from the action Sfi — Sdbi + Swz- The DBI action for the D7-brane is given by 

' 8 



S 



DBI 



j (fx e* 6^ sin e 56 F^ cos' ^ Si S2 S3 , 



(2.10) 
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where we have introduced the following auxiliary dimensionless quantities 

& S3E 



cos^ 



Xq I ■ 2 Xq ^ 

— — -L- cm — — t-t. 



sm 



2 ' F2 2 ' V ^ ' 462 56 ^2 " V ^ ' 52 

with Xq being a function of a determining the brane embedding. The WZ piece of the action is 

S^z^-^Q^ j d^xsm0p{a)b''e^^S^Elcos^^. (2.12) 

The corresponding equation of motion for Xq is given by 



1 + 



1 + -^, (2.11) 



= -d, 



e^cos'^^id^Xq) 
z :i2 



+ 6*6^ S^F^ S3 cos- ^ sm 



2 Xq ■ Xq 

— sm — 
2 2 



+ 



3 .=.1 ^2 



(2.13) 



COS y (^1 - ^ + <3/ e ;^ cos y p(a) ^3 



As we commented in the previous subsection all functions depend only on u, hence it is possible to 
describe the system in terms of a one-dimensional effective action. Inserting all the ingredients in 
(2.6) we obtain: 

,3 



'e// = 2^2 ' '^'^ ' 
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(2.14) 



where Vi^z is the volume of the Minkowski space and Cid is given by the following expression 



1 fh' 



'Id 



2\h 



+ 12 



S' 



+ 



F'S' 



b' fh' ^S' F'\ 1 fh' 
+ i U + '^ + 'f +2 U 



/)2n2 1 

" _ _ n2 «2 1,2 2* r,8 ^2 



(2.15) 



4g/e*6'F'>5^SiS2S3 cos' 



3 



/2 



1 e-* /i 

2 P 



g2$j^ 

H'^ 



QcHJ' . 



Since the potential J enters the effective action only through its derivative, it corresponds to a 
"constant of motion". This new parameter is related to the value of the magnetic field close to the 
boundary through the equations of motion for F3, coming from the lOd supergravity [33]. We will 
fix this constant of motion in the following way 



dCid 



Qc 



(2.16) 



dJ' h 

The next step is to use equation (2.16) to eliminate J', in favor of H^, in equation (2.15) after 
performing the following Legendre transformation 



^id — L 



Id 



5J' 



(2.17) 



J'=J'(H,H^) 
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and then calculate the Euler-Lagrange equations from the transformed action (2.17). Defining the 
following auxiliary (dimensionless) expressions 



^4 = l-COt --, e = COS 

^2 = 1 + — T7?r- & /33 = 1 + 



HI2 



(2.18) 



we can write the equations of motion in the following compact way 



9^(log6) = -AQfR^hS^'F^^ - Q} h /3s , 



(2.19) 



u2 1 /, TT/2 

5^(log/^) = -Qlj^, - ^e'^H'Qjp'hS'Ps + (1-/3^) 



62 



(2.20) 



(2.21) 



^QfeH'F'S' cos^^Sa 



cos2^^ + ^fl-5^ + 2=i 

2 .=1 .=2 V ^3 



(fog F) = 4 6^ - ^3 gj p2 ^2 ^8 ^ 



1 e-'^hH"l3, 
4 P 



(2.22) 



2g^e*62>S« cos^^sin^^^ 
2 2 .=ii.=i3 



1-ls, (1-S^) 



5„ 



62 



1 + 51 



Q) p" e^* ^2 ^8 ^ g<P p-^ ^ 



2 772 



1 

2 P 



e^Q)p^HhS^ + Qc^'+ 4QfHhS^F^^ 



(2.23) 
(2.24) 
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Together with the above system of EOM we get the following "zero-energy" constraint 
2 \h J \ ^ J 2 

+ \q)p' e^^S' Es + AQf 6^ e* ^ cos^ ^ . 

2 ,^^2 2 

Equation (2.25) is obtained by requiring invariance of the one dimensional effective action (2.14) 
under an infinitesimal reparameterization cr — > (1 + 5A)(T. This is equivalent to requiring that the 
one-dimensional Hamiltonian Hid corresponding to Cid vanishes. An important observation is that 
in generating the Hamiltonian Tiid one should Legendre transform with respect to all fields in Cia, 
including the field Xg(cr) specifying the profile of the fiducial embedding. This suggests that we can 
treat x?(c) ^ ^ dynamical variable, therefore we can also obtain the equation of motion (2.13) from 
varying the effective lagrangian (2.15). Remarkably if we use equation (2.5) to substitute for p{a) 
in the effective Lagrangian Cid and derive the corresponding Euler-Lagrange equation for Xg(o'), we 
reproduce exactly the equation of motion (2.13). This provides a non- trivial self-consistency check 
of the construction outlined above. 

The system (2.16) & (2.19)-(2.24) allows for a systematic expansion of all the functions in power 
series oi Qf, as defined in equation (2.4). In fact physically it is more relevant to expand in the 
parameter, e* 

= Q/e*% (2.26) 

which takes into account the running of the effective 't Hooft coupling (through the dilaton factor 
e** in (2.26)). We consider the following first order expansion in e*: 

b=l + e,h, h=^{l + e,hi), S = r{l + e,Si), (2.27) 

F = r (1 + e,Fi) , $ = + e,*i , H = H^{1 + e^Hi) & Xg ^ Xo + e*Xi , 
together with the reparametrization: 

ri^e-^*H',R'' & f=— , (2.28) 

where = Qc/4:. The result is a coupled system of second order differential equations which can 
be decoupled by the transformations: 

Ai = Fi - 5i , Ai = Fi + 45i + ^ 6i & = hi - bi . (2.29) 
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For the decoupled system we obtain: 



Ik- 

r 


^2 A - 


r 




r 


12^ 



1 S2 1 

A'/ + :^a; - ^Ai = - Ja,, a;' + .a; - -Ai = - Ja,, (2.30) 

r r 2 r r I 

where the analytic expressions for the sources , Jh^ , J^i , Jbi , Jai , Jai appear in Appendix A. The 
equation coming from the constraint (2.25) is: 

V'l -Ui- -Uif) (2.31) 
with 771 = 2Ai + Ai and Jr^-^ given in Appendix A. 

As one can see all of the equations of motion, except the one for the fiducial embedding, in (2.30) 
are linear, therefore it is possible to obtain their solution in an integral form in terms of appropriate 
Greens functions. On the other hand the non-linear equation of motion for is the same as the one 
for a probe D7-brane studied in ref. [ours]. In fact it is the only non-linear equation that we have 
in our construction and we solve it numerically. The general solution for the classical D7-brane has 
the following expansion at large f 

«i^^ = ^ + ^ + ^(^^J ' (2-32) 

where rfio = mo/rm & c = c/rj^, while mo & c are proportional to the bare mass and fundamental 
condensate of the dual field theory. Inside the bulk of the geometry xo terminates at a given radial 
distance fmin = at which the ^S"^ cycle wrapped by the D7-brane vanishes. As described at the 
beginning of this section the smearing procedure produces a spherical cavity of radius Vg. Inside 
this cavity the solution is sourced solely by the colour D3-branes through the self dual RR F5 form. 
Before presenting the solution for the functions of the first order perturbative expansion (2.27), we 
will elaborate more on the equations of motion inside the cavity. 

2.3 Vacuum solution inside the cavity 

The smearing of massive flavour D7-branes produces a hollow spherical cavity at the origin of the 
subspace transverse to the colour D3-branes. Inside this cavity the gravitational background is a 
solution of the vacuum equations of motion (setting Q/ = in equations (2.16) & (2.19)-(2.24)). It 
turns out that it is possible to obtain non-perturbative solutions for all background functions, except 
for the functions S and F describing the radii of the CP2 internal subspace and corresponding U{1) 
fiber bundle. 
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2.3.1 The equations of motion in the vacuum 

The equations of motion inside the cavity are given by 

e-^'hH'^ 



62 



dlilogS) = -2 b^F^S^ + 6 b^F^S'' + \ + ^e^ h J'^ , 



9„ 



hH' 



62 



2 62 2 



(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 
(2.39) 



Adding and subtracting in various ways (2.33), (2.35), (2.36) & (2.37) we easily obtain the foUowing 
system of equations without sources 

F2 F2 / p2\ 

-12(65^)^^(1-^), (2.40) 



a^iog 



^2 



52 



52 y ' 



dl\og{bS') = 8 {bS') 



2 F2 / F2 



(2.41) 



92loge-'*6 = . 



(2.42) 



While (2.42) strongly suggests e oc 6 for (2.40) & (2.41) we need to define a new set of variables 



U = bS"" k V = 



52 



In these variables the equations of motion ior V h U decouple in the following way 

dllogU = SW^V {3 -V) , 



(2.43) 
(2.44) 



a^iogy = -i2[/V {1 - V) . 



(2.45) 
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Notice that V is tiie ratio between the radii of the CP2 and the fiber in the S^, therefore it is a 
measure for the relative squashing. The following solution to (2.44) & (2.45) 



V(a) = 1 & U(a) = , (2.46) 

^ ^ ^ ^ 4(7 + const ' ^ ^ 

corresponds to a non-squashed S^. Since in our case the vacuum solution at the boundary of the 
cavity should match the flavour background, we are interested in deformations of (2.46) correspond- 
ing to a squashed S^. We have no reason to expect enhancement of the global symmetry of the 
theory in the deep IR. From a field theory point of view the squashing of the corresponds to a 
breaking of the global SU{4) symmetry down to SU{3) x U{1). Inside the cavity the gravitational 
background corresponds to the effective field theory obtained after integrating out the fiavours. 
Now let us briefly discuss the rest of the vacuum equations of motion. A natural candidate for 
a solution consistent with our ansatz is the supergravity background dual to a non-commutative 
Yang-Mills, which can be obtained as a near horizon limit of the supergravity solution corresponding 
to the D3-D1 bound state. Indeed one can check that upon the following substitution 

^»„^^».l±|ej_ (2.48) 



(2.49) 



j(f) = eV^-^(i + c6e^)-^(i + c$e2)-j, (2.50) 

a{r) = 1^ + c,Q' ; = e^e^Vi?^ , (2.52) 

equations (2.33), (2.34), (2.37), (2.38) & (2.39) are satisfied where ©^^ is the parameter of non- 
commutativity of the dual field theory in the (a;^, x^)-plane. The constants of integration Cb , c$ and 
Ccr will be fixed by matching to the flavour part of the background at the boundary of the cavity. 
Notice that the system of equations (2.44) & (2.45) is completely decoupled and describes the 
different deformations of the compact part of the geometry. Expressions (2.47)-(2.52) together with 
(2.46) constitute a full solution corresponding to the gravity dual of a non-commutative SYM with 
SU (4) global symmetry. It is possible to flnd a more general class of solutions corresponding to a 
non-commutative SYM with SU{3) x C/(l) global symmetry and this will be the topic of the next 
subsection. 
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2.3.2 Deforming the sphere 

We now construct a non-commutative non-supersymmetric background by considering perturbative 
solutions of the system (2.44) & (2.45) around (2.46) starting from the following ansatz: 

t/((j) = ^ (1 + Emi) & V{a) = l + Evi, (2.53) 

where ui & vi arc functions of a to be determined after substituting in (2.44) & (2.45) and expanding 
in E. At first order we have 

< - ^ ai- ^ & v'l-^Vi^O, (2.54) 



a 



and after an appropriate redefinition 

P=^ui + vi ^ /3"-^/3 = & v';-^vi = Q. (2.55) 

The general solution is given by 

P{a) = cicr^ + C2(7~^ k vi ^ Csa^^'^ + c^a'^l'^ , (2.56) 

and in terms of the radial variable r defined in (2.52) the functions ui & vi are given by 

1 2 2 

Ul = Ci — + C2f^ - -csf"^ - -C4f^ , (2.57) 
r'^ 5 5 

^^1 = cgf-^ + c^P , (2.58) 

where the constants c, are related to the constants Cj by a non-singular diagonal linear transforma- 
tion. Notice that in the formalism of (2.53) the parameter of non-commutativity is proportional 
to E and hence a = l/(4f*^) + 0(E). Analyzing (2.58), the general solution regular at f = is 
spanned by C2 7^ 0, C4 7^ , ci = C3 = 0. The constants C2 and C4 will be determined by matching to 
the flavour solution at the boundary of the cavity (after setting E ~ e*). 

2.4 Constructing the full perturbative solution 

Our construction is a perturbative expansion near a supersymmetric AdSs x background. Re- 
quiring that the leading order corrections inside and outside the cavity agree implies that ©■^ ~ e*. 
Furthermore, we introduce the radial variable r related to a via = ^ and its dimensionless analog 
f defined in equation (2.28). The resulting expansions are 

(2.59) 
(2.60) 

(2.61) 

(2.62) 

(2.63) 



b = 


1 + CbO" , 


$ = 


+ - c<^ e^ 


H{r) = 
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Since (2.61) docs not depend on any constants of integration we can use the perturbative expansion 
of the flavour solution to determine the exact relation between and e*. Indeed comparing (2.61) 
and (2.27) we conclude that for f < fg we have 



2„% ~4 



i?2 



(2.64) 



where is the radius of the cavity. In the following subsections we will present semi-analytic 
solutions for all the functions of the first order perturbative expansion, requiring specific behaviors 
inside the cavity emerging from the expansions (2.59)-(2.63). 

2.4.1 Solving for the Meld 

In order to solve the equation of motion for Hi in (2.30) we need the Greens function G'i?i(r, fi) 
satisfying the following equation 

-^G(f,fi) + ^-^G'(f,fi) - ^G(f,fi) = 5(f-fi) 
ar^ r dr 



(2.65) 



The solution for Hi that inside the cavity is ex and vanishes at f = f* can be written in the 
following integral form 



Hi 



dfi JHiih) 



GHi{f,fi) - ^GHi{f^,fi) 



(2.66) 



Elaborating on (2.66) inside the cavity (f < fq) we have 

Hi{f)^-H,^ (f„f*) f^cxf\ 
where the expression for the constant # {fg, f*) is 



(2.67) 



(2.68) 



Comparing (2.64) and (2.67) we obtain the promised relation connecting the non-commutativity 
parameter to the number of flavours 



Nf 1 



(2.69) 
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2.4.2 Solving for the dilaton &; bi 

In order to solve the equations of motion for both the dilaton and 61 in (2.30) we need the Greens 
function G'{<i.i,bi}(f, fi) satisfying the following equation 

^G(f,fi) + ~^G'(f,fi) = S{f-h) 



G'{$i,6i}(r,fi) 



I (1- Jl ) 



(2.70) 



Looking at (2.59) and (2.60) wc sec that both the dilaton and b are constant inside the cavity. The 
solution for $1 & bi which is constant inside the cavity and vanishing at f = can be written in 
the following integral form 



G'{*i,6i}(r,ri) - G'{$,,6i}(f,,fi) 



Mr) = 2 / dfiJbAh) 



(2.71) 



(2.72) 



(2.73) 



(2.74) 



Inside the cavity both $1 and hi are constant and are given by the following expressions 

ft 

^i(r) = ~\ J ^^1-^*1(^1) G'{*i,6i}(^*>^i)> 

ft 

bi{f) ^ - \ J dfiJbAh)G{^,^b^y{h,fi) . 

Equations (2.73), (2.74) and (2.69) can be used to fix the constants Cb and c#. 
2.4.3 Solving for Ai 

Proceeding in the same way we obtain the Green's function for the equation of motion for Ai, defined 
in (2.30) : 

-^G(f,fi) + ^-^G(f,fi) - ?^G(f,fi) = 5{f-fi) 
dr'^ r dr 

~4 



Gx.{r,h) = (1 - ^ ) - ' (2-75) 
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and look for a solution that inside the cavity is oc and vanishes at f = f* 



(2.76) 



2.4.4 Solving for Ai 

Proceeding in the same way for Ai we have 

-^G(f,fi) + ^-^G'(f,fi) - ^G{r,ri) = 5{r-ri) 
ar^ r dr 



GA,(f, fO = ^ (l - ^) e(f-fO. (2.77) 



Looking for a solution which behaves as Ai(f) oc P inside the cavity and asymptotes to the 
supersymmetric solution at f* 



(B.13) 



(2.78) 



we have 



Ai(f) = - / dhJ^,{h) 



^2 



GAi(r,fi) - T^GAi(r*,ri 



12 f2 



1 - 



2\ 3 



(2.79) 



2.4.5 Solving for Ai 

The Greens function for Ai is given by (2.75) so we proceed immediately in looking for a solution 
which behaves as Ai(f) oc inside the cavity and asymptotes to the supersymmetric solution at 



(B.14) 



Ai(r* 



m 

' 18 I " ^.2 



3 



we have 



Ai(0 = 1 J dhJA,{h) 



~4 



GAi(f,fi) - T7G'Ai(r*,ri) 



1 f-^ , 



;^2\ 3 



m, 



(2.80) 



(2.81) 



2.4.6 Solving the constraint for % 

The first order differential equation for r} (2.31) can be integrated to give 



+ dr — ^- — 9 r - Tq 



with 771 (r*) =^ ^ I 1 — 



;^,2\ 3 



(2.82) 



The value of r]i{f*) is fixed by the requirement of obtaining the supersymmetric solution at f = f*. 



18 




Figure 2: Family of solutions of the background functions for different bare mass parameter m. 



2.5 Numerical solution 

As we learned from the previous subsection it is possible to obtain a solution for our supergravity 
background in a semi-analytic form in which all background functions are written in an integral 
form. The solution is completely determined by the D7-brane charge distribution function p{r) or 
equivalently by the profile of the fiducial embedding Furthermore, in the perturbative expan- 



19 



sion the equations of motion for the first order corrections to the background fields are sourced by 
the zeroth order expansion of x?(^) (namely Xo(^)), which is known from the probe approximation. 
The equation of motion for Xo(^) is non-linear and we solve it numerically. 

An obvious approach to obtain a numerical solution for the background functions would be to solve 
numerically the integral expressions presented in section 2. Actually it is more convenient to solve 
the equations of motion completely numerically using Mathematica's built in function NDSolve. To 
generate the solution we employ standard shooting techniques and sew the numerical solution to 
the analytic solution inside the cavity at f = f g . The constants of integration (which specify the 
shooting parameters) are obtained by matching the solution to the supersymmetric one at f = f^, . 
This approach has the disadvantage that the accumulated numerical error grows strongly with the 
parameter f*, which we need to keep numerically high. It can be improved if we approximate the 
function xo with a taylor series expansion in 1/f for f > fsew {fsew is a number of order ten) and 
solve analytically for the background fields in the region fgew < ^ < ^* by matching the solution to 
the numerical one at fsew and to the supersymmetric one at f^,. 

Repeating the procedure outlined above for different values of the radius of the cavity fq generates 
solutions corresponding to different quark masses. In figure 2 we have presented plots of the back- 
ground functions for the following range of the bare mass parameter mo £ [0, 2.3], which is the one 
needed for the analysis in section 5. As one can see the solutions corresponding to different values 
of fh differ significantly in the infrared f ^ and are very close to each other for r < f*. In fact 
one can check that, except for the function Hi, for f < f^, the background functions are very well 
approximated by their corresponding functions from the supersymmetric limit. The function Hi 
follows the pattern ifi(f) ^ ^(?^ — 1) represented by the black dashed curve in figure 2. Appar- 
ently at f*, i?i(f*) vanishes, but has a non- vanishing first derivative, therefore the matching to the 
supersymmetric solution is not smooth at f*. However, one can check that i?((r*) = 1/f* which is 
sub-leading for f^^l. 

Overall we conclude that our solution is very well approximated by the supersymmetric background 
near f* and hence one can approximate the background with the supersymmetric one for f > f*. 
This is particularly useful to investigate the UV behaviour of the background, because it relates 

non-pertTirbativcly the arbitrary UV scale and the parameter of pcrturbativc expansion to the 
position of the Landau pole Tlp. In the next subsection we provide a detailed description of the 
hierarchy of scales and regime of validity of our perturbative solution. 

2.6 Hiercirchy of scales and regime of validity 

In this subsection we analyze the regime of validity of our perturbative solution and the hierarchy 
of energy scales (in terms of radial scales) of the theory. Our analysis follows closely section 2.4 of 
ref. [29], where the finite temperature system has been unquenched. In terms of radial coordinates 
the hierarchy of scales can be written in the following way: 

< ~ < r* < r„ < tlp , (2.83) 

where and Vq represent IR energy scales related to the energy scale set by the magnetic field 
and the one corresponding to the physical mass of the fundamental matter. The radial scales Va 
and tlp represent UV energy scales corresponding to the scale at which the supergravity solution 



20 



develops pathologies and the Landau pole of the theory at which the effective 't Hooft coupling 
blows up. For energy scales close to r* our solution, as can be seen from the plots in figure 2, is well 
approximated by the supersymmetric solution corresponding to vanishing magnetic field [28] , which 
(given our choice <C r^,) is well approximated by the supersymmetric solution corresponding to 
massless fundamental fields [27]. This is why the analysis of the UV energy scales is exactly the 
same as in ref. [29] and in particular the relation between the "position" of the Landau pole tlp, 
the finite cutoff r* and the perturbative parameter e* is: 

— e-^l^* < 1 . (2.84) 

Furthermore in order for our perturbative solution to be valid in the region < r < we 
need e*(''«''/e** ^ 1, which requires |$i(fg)| <^ 1 (because |$i(fg)| is the biggest off all first order 
corrections). However, as one can check from the plot in Figure 2, maxj=^|$i(fg)| ~ log y. Therefore 

we need ^ > e~^/^*. 

In addition we would like to be able to neglect terms sub-leading in the r^/r* expansion. Therefore 
we have to make sure that the corrections e^, that we are considering are much larger than the one 
that we ignore. This requires e ^ ^. In summary we have: 

Auv n Sn^Nc 

Finally, validity of the supergravity approximation requires that we ignore closed string loops (A^^^ 3> 
1) and a' corrections (A^ ^ 1), where Xq is the effective 't Hooft coupling at the energy scale set by 
Tq ~ Mq. It is related to A* via Xq — A». In addition, validity of the smearing approximation 
suggests large Nf. In summary we have: 

7V,»1, A,»l, e,= ^^«l, (2.86) 

where we have defined the IR perturbative parameter Cq. Clearly it is related to (defined in 
equation (2.26)) via: 

= = e*(l + e,$i(fg)) + 0{el) , (2.87) 

e * 

which also implies: 

e, = e,(l-e,$i(rg) + 0(ej) . (2.88) 

Finally requiring that a' corrections, which scale as Ag^*^^, (for more details look at ref. [29]) are 
sub-leading relative to flavour corrections controlled by requires: 

A-'/' < eq . (2.89) 

We close this section by commenting that the numerical values used in our analysis are in the regime 
of validity specified above. 
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3. Free Energy and Condensate 



In this section we calculate the free energy density and fundamental condensate of the dual field 
theory directly from the supergravity background. 

3.1 Helmholtz versus Gibbs free energy 

Following the general prescription of [39] we identify the on-shell Euclidean action with the Helmholtz 
free energy. The Euclidean action has contributions from two terms Xbuik and Xgurf given by: 

^buik = ~ 2 [ ^UBda , (3.1) 

I.uri = -^VlK , (3.2) 
'^10 

where Tub is the wick rotated action (2.6) and in (3.2) is the standard Gibbons-Hawking term. 
In fact one can check that the one dimensional effective lagrangian £id defined in (2.15) already 
includes this boundary term. Therefore we can write: 



^hulk ~l~ -^surf 



^ / Cu{B,B',J',xl;,Mda , (3.3) 



where ipa is a collective notation for the background functions (6, h, F, S, $). Note that: 

ST dCid , . 

Jjcc — CH., (3.4) 

where we have used (2.16). Given that we make the identification X/V4 = J^/V^, where J-' refers 
to the Helmholtz free energy and we have that |^ = i?*, it is natural to relate the field J to the 
magnetization M of the system. However we are interested in a thermodynamic ensemble in which 
we keep the external magnetic field if* fixed. Therefore the proper thermodynamic potential is 
given by the Gibbs free energy Q = — MH^. Equation (3.4) implies: 



X= - "^J Lu{B, B', H,,,jja, i^'a) da , (3.5) 



were Cid refers to the Legendre transform of Cid defined in (2.17). Next we expand in e* using 
(2.27) and define dimensionless variables along (2.28). To first order one has: 

X = To + e^XoBi + eAound + 0{el) , (3.6) 



and more explicitly: 

'-g = 6rt - 4e*r^ /^^ " J? (r^'vir)) + rVl + f^i{r) 



9k2 



(3.7) 



6r^ + 4e*r^77(r*) - 4e* J drf Vl + r^^{r) + 0{el) 
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The first two terms in (3.7) can be cancelled by appropriate counter terms (we refer the reader to 
the next section for more details) , while the last term is precisely the DBI-term describing a probe 
D7-brane. Therefore to first order in ~ Nf/N^ the Gibbs free energy of the unquenched system 
coincides with the free energy calculated in the quenched approximation. This feature has been 
observed also in the finite temperature case studied in [29]. This suggests that to first order in 
one cannot study the effect of the dynamically generated energy scale set by the Landau pole of the 
dual gauge theory. Therefore we need to compute the second order contribution to the Gibbs free 
energy^. Furthermore, despite the fact that we arc dealing with a finite cut off set by the parameter 
r^, we can still perform a holographic renormalization of the free energy resulting in an expression 
finite in the large limit. 



3.2 Holographic renormalization 

In this subsection we will regularize the Gibbs free energy defined in (3.5). One can show, that 
modulo logarithmic divergences, the counter terms needed to regulate (3.5) are the same as ones 
needed to regulate the supersymmetric background corresponding to vanishing magnetic field (the 
— limit). Therefore, it is natural to regulate the on-shcll action in (3.7) by subtracting the 
supersymmetric on-shell action. Note that this would suggest regulating the DBl contribution to 
the free energy by subtracting the DBI term for the supersymmetric case. Such an approach would 
make comparison to the results obtained in the quenched approximation difficult, where the on-shell 
action is regulated by the addition of the appropriate boundary terms. Furthermore, logarithmic 
divergences depending on the magnetic field do not have analogues in the supersymmetric back- 
ground. To rectify this we consider a mixed regularization scheme: we regulate logarithmic and 
divergences due to the DBI term by adding appropriate covariant boundary counter terms, while 
the rest we directly cancel by subtracting the relevant part of the supersymmetric action. 



3.2.1 The supersymmetric case 

In the supersymmetric limit of vanishing magnetic field, equations (2.13) and (2.19)-(2.24) can 
be integrated to give a BPS system of first order differential equations^. The full non-perturbative 
solution for the corresponding supersymmetric background was obtained in [28] and for more details 
we refer the reader to the appendix B. Here we present the expansion to second order in e* of the 



^Note that from a holographic point of view even at first order in e« our ansatz provides a novel feature, namely 
the relation between the field J(r) and the magnetization of the dual gauge theory. 
®Note that in this hmit 6=1 and H = J = 0. 
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supersymmetric Euclidean "on-shell action" and analyze its divergences. We obtain: 



2 k 



10 . 



3-^susy 



6r! + 4e* 



drr^ cos^ Y\^^ 



mo 



+ 



(3.8) 



cos^ ^X'o 



r=r* 5 



mo 



mo 



mo 



+ 0{el) . 



As we pointed out above the divergences of the action (3.8) completely cancel (modulo logarithmic 
terms) those of the non-supersymmetric action. However we would like to regulate together the 
DBI term in the first order expansion and the term depending on Xi in the second order expansion 
(coming from expanding X = Xo + e*Xi i^i the DBI term) by adding covariant counter-terms. We 
make this choice in order for our regularization scheme to be compatible with the one employed 
in [13], for the probe approximation. Furthermore, by construction, boundary terms which do not 
contain derivatives of the fields are exactly the same with those of the non-supersymmetric action 
and can be directly subtracted. We define a subtracting action: 



f — '7-bdr I 'T-bulk 
-i-subt — -^subt "T -^subt ' 



(3.9) 



where: 



^'^lO T-bulk _ 2 



3-''subt 



+ 



J drr' {Jl^Al - J^.Al) 



mo 



(3.10) 



^ f drr^ cos' Y ~ ^ ^ drrV|$^+ 



mo 



mo 



162^^! 



(3.11) 
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and we have explicitly evaluated I^^^- Next we proceed with the regularization of: 



%usy + X,.bt) = 46. fdrr' cos^ ^ \^^^^ + ef-^^^£Mxi + 0{el) 

J Z \ ^ /-, , /2 r=r* 

mo Y 4 AO 



3 



4e*y" drr^ cosily 1 + ^x" + C(e*)- (3-12) 



mo 



Note that in terms of the full function x(r) the subtracted supersymmetric action looks just like 
the DBl term of the probe approximation. However, the function x(r) is not the same (namely 
2 arcsin(m/r)). 

The K-symmetry condition for the embedding can be solved in a new radial coordinate f(r): 

f{r)^r{l + e,pi{r) + 0{el)) . (3.13) 

where pi(r) is given in (B.7), with pi(r*) = and hence f(r*) = r*. It is in this new radial coordinate 
that the fiducial embedding satisfies: 

x{r) = 2arcsin— , (3-14) 

where m — f(mo). In fact we identify m as the full bare mass parameter m = f(mo) — mo + e*mi + 
0{el). Clearly this defines mi in terms of mo- Furthermore, since f(r*) = r* one can write: 

x(n) = ^ + ^ (^) = xoin) + e.xi(n) + O {el) = ^ + e^xiin) + ^ (^) + . (3-15) 
suggesting the natural relation between mi and the leading term in the expansion of Xi(^*)" 

X,W = ^+o(^). (3.16) 

Equation (3.16) is an alternative definition of the first correction mi to the bare mass, which does 
not require the existence of a special radial variable f. We will use this definition when we study 
the theory at finite magnetic field. 

The fact that at r* the fiducial embedding still has the same expansion as in the quenched approx- 
imation suggests that one can use the same counter terms. Indeed one can check that the following 
counter terms: 



(3.17) 
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completely cancel the first order term in (3.12) and leaves only sub-leading terms at second order 
in e*. More precisely: 

-^susy ~l~ -^subt ~l~ -^^ount ~l~ ^* 

If the theory were UV complete the addition of the regulating terms Xgubt a-nd Xcount to the "on-shell" 
action would provide a vanishing mass dependent expression for the free energy of the supersymmet- 
ric background. However, due to the existence of a Landau pole our UV cutoff is finite. In general 
this can lead to spurious mass dependence of the free energy (and hence spurious condensate) due 
to the #(m) term in (3.18). Fortunately a similar contribution from sub-leading terms is present in 
the non-supersymmetric (finite magnetic field) case too. One can verify that the two contributions 
are approximately equal for large bare masses, which is also the regime when these contributions 
are significant. Therefore, ignoring sub-leading terms is a very good approximation, because they 
can always be canceled by a redefinition of Xgubt- In the next subsection we will regularize the free 
energy at finite magnetic field by adding IsuU + Xcount to the action I. 



+ O(#(m)^ 



+ O {el) 



(3.18) 



3.2.2 Regulcirization at finite magnetic field 

In this subsection we regularize the "on-shell" action I by adding the term Igubt defined in (3.9) and 

the counter terms Xcount defined in (3.17). It is convenient to split the contributions to I into two 
parts I = + Xbuik, regulated by Xcount and Xgubt respectively. We find: 



2 k 



10 



V4TT^ 



Xbulk + Xsubt) — e*r: 



2^4 
m 



114 1 
^Jbih + ^JxAi - g (^AiAi - JaiAi) + - 



+ 



(3.19) 



drr 



Jh.Hi +{f^ + ^) cos' 



Xo 



--rt(AiA/-AiA/) 
o 



rl-miy{nrl + 7ml) 
162r! 



The expression in equation (3.19) can be evaluated numerically and analytically in the limit <^ 
mo (weak magnetic field). In the limit of weak magnetic field one finds that the subtracting term 
regulates all divergences. We will study this in more details in the next subsection. Let us now 
focus on the remaining part of the action: 



9k2 



{■^X ~^ -^count) ^* 



I 

4r^ j drrVl + f^ cos^yyi + 



-Xo 



/2 _ ^4 / 2 



Xo 5xo 
2 48 



+ 



f^VT+f^cos^ ^x'o 



1 + tXo' 



-Xo + 



5Xo 
12 



Xi{r*) ■ 



(3.20) 
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The first term in the first order contribution to the action + Xcount is the same as in the quenched 
approximation and hence one needs to complete the counter term action to subtract the oc log r* 
divergence. The complete counter term action is: 



(-^count JDaBJD 



"10 -7- 



Now using the fact that Xo(^*) = 2mo/f* + 2co/fJ + 0(l/f^) one obtains: 

' log r 



log 



R 



(-^x -^count) — 



^D7{rho) + e* -2co mi + O 



+ 0{el) , 



(3.21) 



(3.22) 



where we have defined 2D7(m) and rhi via the following relations: 



^07(^0) = y c? f fv 1 + f'^if (f) 



/- X 2mi /^logr, 
XU?^*) = -^ + 



~4 



4 V 2 ^ 48 J- 2^ R' 



(3.23) 
(3.24) 



The function 2D7(mo) is proportional to the expression for the free energy of the fundamental matter 
in the quenched approximation. It is known that (see for example [40]) dLx^-jimo) j d-m^ — — 2co. 
Using this relation (3.22) can be written as: 



9k2 



(-^X -^count) — T, 



lB7{rho) +XD7(mo)(e*mi) + ( O 
logf, 



log 



+ 0(e^) = 



4e,r^XD7(m) + ef ( O 



+ 0{ei) , 



(3.25) 



where we have defined the complete bare mass parameter as m = mo + mi + O(e^). Note that 
modulo sub-leading terms the only contribution from the el term in (3.22) is the one needed to 
complete the argument of Xd7. Now we can write our final expression for the regularized action 



-^reg — X -|- Xgy^j^ -|- X(.Qmit ■ 



^^reg - 4e. rt i^rirh) + el r^xg (m) + 0(6^) , 



(3.26) 



where the quantity X^j[m) is given by the right-hand side of (3.19) divided by e^r^ and we have 



completed its argument to a full bare mass m. We are now ready to proceed with the evaluation of 
the free energy and the fundamental condensate as a function of the bare mass. We will study first 
the limit of weak magnetic field (m ^ 1), which can be analyzed analytically. 
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3.3 Free energy and condensate at weak magnetic field 

In this section we calculate analytically the Gibbs free energy Q and the fundamental condensate of 
the theory in the limit of weak magnetic field. Given that if* oc and rh — m/rm this corresponds 
to the 1/m expansion in equation (3.26). It has been shown in [13] that: 

XD7(m) = logm + consti + 0{l/m^) . (3.27) 

Furthermore, using the approximate analytic solution for the background (see Appendix C) one can 
show that: 

X^f)(m)= const. + (3 

Note that, since we arc dealing with a finite cutoff, numerically, the logarithmic term in (3.28) is of 
order one and hence the sub-leading term is of order All together we obtain: 

|^X,eg(m) = -2e* (r^ log m + const) + 0(e2) oc -e42//,Mogm + #(//*)]+ 0(e2) , (3.29) 

where we have ignored terms of order (1/m*^). For the leading order contribution to the free energy 
we recover the result from the quenched approximation: 



1 . TT^ 

12" 



^-g{H,,m) = -e,^ri\ogm + 0{el) , (3.30) 

V3 



where we have skipped the mass independent term #(if*). Now we can define the fundamental 
condensate via: 



, ISSIH^ = + 0(H;) + . (3.31) 

V3 dm kJq m 

Equation (3.31) suggests that to order the fundamental condensate is the same as in the quenched 
approximation. However, as commented in section 4.6, the relevant perturbative parameter which 
takes into account the running of the effective coupling is e^. Using (2.88) we arrive at the following 
result for the condensate: 

^(^^) = -e, ^ (1 + + log ^)) + 0(4) , (3.32) 

where we have used that for fg ^ 1 to leading order one has: 

$i(f,)«.-^-log^«.-^-log^. (3.33) 

Therefore our conclusion is that at first order in the result is the same as in the quenched 
approximation. However, at second order in tg we observe a logarithmic running of the condensate 
as a function of the bare mass and the energy scale set by the finite UV cutoff Auv ~ r*, which 
reflects the positive beta function of the theory. We proceed with calculating the Gibbs free energy 
and the fundamental condensate at strong magnetic field. 
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3.4 Free energy and condensate at strong magnetic field 

In this subsection we calculate the Gibbs free energy and the condensate of the theory at strong 
magnetic field evaluating numerically the functions XDi{m) and Zf^j{m). As commented above the 
first function describes the Gibbs free energy to first order in and has been studied in [11, 13]. 
Apparently it also describes the first order contribution to the fundamental condensate given by: 



10 



{^P^|J) = e,i^,,{m) + 0{el 



-2e, do{m)+0{el) 



(3.34) 



In figure 3 we present plots of Xj^Tim) and — Co(m). As one can see from the plot of Xj^jim) for 




m 



Figure 3: Plots of Xd7 and cq versus the bare mass parameter rh. 

sufficiently small fh the theory has multiple phases. However only the lowest branch is thermody- 
namically stable and studies of the meson spectrum in the quenched approximation [9] verified this 
stability. In the second plot of Figure 3 we see the condensate, but only for the stable phase. At van- 
ishing bare mass the theory has a negative condensate, while for large bare masses the condensate 
has the — l/(4m) oc Hl/m dependence governing the weak magnetic field limit. 

Our next task is to obtain similar plots for the second order contribution in e^, to the Gibbs free 
energy and the fundamental condensate. We will restrict ourselves in studying only the stable (to 
first order in e^) phase of the theory. Using the numerical solution for the background we study 
X^^;(m) and its derivative. Note that (3.26) and the definition of the condensate: 



1 dQ{H^,m) _ 1 dXrcg{H^,m) 



dm 



dm 



(3.35) 



imply the relation: 



10 



27r^ri' 



(ijij) = tS'uiim) + ^tlx'^]'{m) 



+ Oie: 



-2e* Co(m) 



2e%{m) 



0{el) 



(3.36) 
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L5_ 2S) 




~ (2) 

Figure 4: Plots of Ify^j and — C2 versus the bare mass parameter m. 



where C2 is defined by Xp^' (m) = — 8c2(m). The corresponding plots are presented in figure 4 
and from these one can see that C2 approaches zero much faster than cq in Figure 3. Surprisingly 
the second order contribution to the fundamental condensate is positive. However the relevant 
perturbative parameter is (defined in section 2.6) therefore in order to obtain the second order 
correction to the condensate we need to trade 
^i(^g) < (see Figure 2) we obtain: 



for eg in (3.36). Using (2.88) and the fact that 



271-3 

m 



{^^) = -2egCo(m) -2e^(£2(m) + \^{f,)\~co{m)) + 0{e^ 
= -2egCo(m)-2e2(c2(m)) + 0(e^) , 



(3.37) 



where we have defined: 



C2[m) = C2[m) 



|$(rg)|co(m) 



(3.38) 



As one can see from figure 2, |<I>(fg)| > 1. Furthermore, comparing Figure 3 and Figure 4 one can 
see that |co(m)| > |c2(rfi)|. Therefore we conclude that — C2(m) < and the second order correction 
to the condensate is negative. This implies that taking into account internal fundamental loops in 
the calculation of the condensate does not change qualitatively the results obtained in the quenched 
approximation. In fact it is not difficult to obtain a plot of C2(m). Our results are presented in 
Figure 5. One can see that indeed the second order contribution to the fundamental condensate 
is negative. The black dashed curve represents the function -^{^ + log^), which governs the 
asymptotic behavior of the condensate at large fh (weak magnetic field). 

Overall our conclusion is that to leading order in a perturbative expansion in the ratio Nj/Nc, 
the free energy and fundamental condensate of the theory agree with the results obtained in the 
quenched approximation. Furthermore, at second order in Nf/Nc we observe that the effect of 
magnetic catalysis is enhanced and the contribution to the condensate of the theory from internal 
fundamental loops runs logarithmically with the finite cutoff Ajjv 
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Figure 5: A plot of — C2 versus fa. 



4. Summary &; Conclusions 

Let us briefly summarize our conclusions and outline possible directions for future studies: 

In section 2 of this paper we have constructed a backreacted supergravity background holographi- 
cally dual to an SU(A^c) A/" = 4 SYM theory coupled to Nj flavours of A/" = 2 hypermultiplets in the 
presence of external magnetic field H^,. Our solution is perturbative in a parameter which counts the 
number of internal fundamental loops. At small distances the geometry has a hollow cavity, where 
it is similar to the supergravity dual of non-commutative SYM theory. The radius of this cavity, r^, 
sets the energy scale corresponding to the physical mass of the fundamental fields. From holographic 
point of view, the supergravity solution inside the cavity corresponds to the low energy effective 
field theory obtained after integrating out the massive flavour fields. The non- commutative nature 
of the theory reflects that fact that the fundamental fields are coupled to an external magnetic field, 
H^,. We have shown that the parameter of non-commutative along the plane perpendicular to the 
magnetic field, G^^, scales as G^'^ 



At large radial distances, our solution is well-approximated by the non-perturbative supersym- 
metric solution corresponding to a vanishing magnetic field [28]. We introduce a sufficiently large 
radial parameter (r* ^ r^) at which we match our perturbative solution to the supersymmetric 
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one. In the holographically dual field theory, the radial parameter corresponds to a finite cutofi^ 
Auv At a given radial distance r^p (r^p ^ r* ), the dilaton field diverges, refiecting the positive 
5-function of the dual gauge theory and the existence of a Landau pole at an energy scale, A^p 
(A^p ~ Tip). We have the following relation [28], oc log implying that to keep the finite 

cutoff sufficiently far away from the Landau pole we need to keep the parameter e* ~ ^ sufficiently 
low. Clearly in the limit — > wc can make the cutoff A^/y arbitrarily large, on field theory side 
this corresponds to the quenched approximation. In fact our studies of the free energy and the 
fundamental condensate of the dual gauge theory confirm that at leading order in e*, we reproduce 
the results of the quenched approximation. 

In section 3 we have applied our construction to the phenomenon of magnetic catalysis, i.e. of mass 
generation in an external magnetic field, beyond the quenched approximation. Note that our studies 
arc at fixed value of the external magnetic field, H^. This implies that the relevant thermodynamic 
potential is the Gibbs free energy of the dual gauge theory, as opposed to the Helmholtz free energy. 
In our holographic setup, the Gibbs free energy is given by the Legendre transform of the "on- 
shell" action with respect to the magnitude of the Ramond-Ramond two-form at r*, J{r^). To 
calculate the free energy, we have developed an appropriate holographic renormalisation scheme. 
This scheme is designed to facilitate the comparison of our results to the results of ref. [8, 11, 13] 
obtained in the quenched approximation. To this end we regularize the contribution to the free 
energy which due to the DBI term of the smeared fiavour branes by the addition of appropriate 
covariant counterterms. The remaining contributions to the "on-shell" action are regularized by 
subtraction, using the supersymmetric solution corresponding to the limit of a vanishing external 
magnetic field [28] as a reference background. 

The application of our regularization scheme to the supersymmetric background [28] implies 
the identification of the bare mass parameter of the fundamental fields m with the leading order 
coefficient in the large expansion of the profile of the fiducial embedding x(r), namely x(r*) = 
2m/r^ + . ... In the limit — )> {r^ — )> oo), this definition of the bare mass parameter m agrees with 
the definition in the quenched approximation given in ref. [24] . Moreover we define the fundamental 
condensate of the theory as the derivative of the Gibbs free energy with respect to the bare mass. 
At first order in e*, our results agree with the results obtained in the quenched approximation. 
At next order in Nf/Nc, the relevant perturbative parameter is the IR parameter defined via 
eg = e'^^^'^y e'^*€^ . Note that takes into account the running of the effective 't Hooft coupling, 
encoded in the running of the dilaton field^. Our results for the second order in Sg contribution to 
the fundamental condensate show that the effect of magnetic catalysis is enhanced. Furthermore 
we observe a logarithmic running of the condensate as a function of the UV cutoff of the theory, 
reflecting the positive /3-function of the theory. This is the main result of our study of the effect of 
magnetic catalysis beyond the quenched approximation. 

An obvious direction for future studies is to investigate the effect of finite temperature. From 
supergravity point of view, this corresponds to generalizing our perturbative construction by sub- 
stituting the AdSs x zeroth order supersymmetric background with the AdSs x black hole. 

^Note also that to first order in e*, we have Cg = e* -|- O (e^) . 
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Such a study would enable us to study the phase diagram of the theory, which in the quenched 
approximation was analyzed in refs. [11, 13], beyond the quenched approximation. It would also 
be interesting to generalize our holographic setup by turning on a finite chemical potential and 
to investigate the effect of internal fundamental loops to the quantum critical points reported in 
refs. [18, 19]. 



Finally it would be interesting to generalize our holographic setup to other Dp/Dq-brane in- 
tersections and to apply our framework to related phenomena also taking place at finite magnetic 
field. For example, it would be interesting to construct the analogue of the effect of Inverse Mag- 
netic Catalysis [41], observed in the Sakai-Sugimoto Model, for holographic gauge theories dual to 
the Dp/Dq-intersections and to take into account the backreaction of the flavour branes. Other 
effects of potential interest include the Chiral Magnetic Effect [42] , investigated holographically in 
refs. [43, 44], and the recently proposed effect of superconductivity from p meson condensation in 
the QCD vacuum in a strong magnetic field [45], investigated holographically in refs. [46, 47]. 
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A. Analytic expressions for the calculation of the EOM and the Free 
Energy 

The analytic expressions for the sources appearing in (2.30) and (2.31) are 

4|(r) , 8 1(f) 



fVf^ + 1 ' 



, 1 



fVf^ + 1 ' 



8f(f) 



JA,(f) = 1(f) 



\/rT 



^4 



1 



3 / 5 . 

-? — - 1 — - COSYnlr, 



with 



e(f)^cos^^yi + -^,{ff. 



(A.l) 



(A.2) 



B. Useful ingredients from the massive supersymmetric solution 

The massive supersymmetric solution was found in [28, 29] and it is easy to check that the following 
set of first order differential equations 



3 _ 2 — - ^ cos^ ^ 



(B.l) 



Xq 



■'4 „^„4 Xq 



d^S ^ S^F\ d„Xq = -2 tan ^ , 9,$ = e* cos^ ^ ■ 

together with 6 = 1 & i7 = solve the full set of equations of motion, (2.13) & (2.19) - (2.24). 
Changing coordinates from a to p, according to dp — da, we can integrate immediately the 
equation of motion for the embedding 

V pPl 
A ?/)7J ^ 

~ eP 



sm 



(B.2) 



where pg is an integration constant related to the bare quark mass. Substituting (B.2) in (B.l) we 
obtain the following set of explicit expressions for 5, F, $ & /i for p > pg 



S ^ e'[l + e^As]-^ , 

$ = - log [1 + e* A$] , 



F = 



dph 



[1 + e^A^]^ 
' [1 + e,As]'^ 



(B.3) 
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with 



1 
4 



^4pg-4p, 



(B.4) 



6 6 2 4 4 



where is the UV scale where the dilaton takes the value and = Qf e** . Following the 
analysis of the main part of the paper, we redefine the radial variable in such a way that the warp 
factor keeps the standard AdS form 



h 



with = -Q^ , 



(B.5) 



and expand dph from (B.3) in powers of e*. Integrating, we obtain an expression for r as a function 
of p, while we fix the integration constant by requiring r(p*) = r* = e''*. Inverting that relation 
we have 



p = log r + e* pi , 



(B.6) 



with 



Pi 



1 

720 



1201og- + 101-- 1-3-^ + 



rrin 



1 - 



ao 



,10 



(B.7) 



+ 120 



rrin 



15^ 



Substituting (B.6) and (B.7) in (B.3) we have 



S ^ r{l + €,Si) , F = r (1 + Fi) , 



$ = + & Xo = 2 arcsin — + e*Xi , 

r 



fB. 
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with 



Xi 



1 

24 



1 + 



1 



1 

24 



1 - 



= log - + 



3r4 

1 

3^ 



+ 



1 TTiQ 

6 



3 rriQ 
16 



1 - 



1 

9 r! 



15 



1 - 



1 r 
6 r 



10 



10 



J^m| lm^/ lr8 
12 r2 16 H V 3 r! 



10 



60 r6 V 3 



4H 



mo 



360 - m2 

+ 120 



- 120 In 



mo 



+ 81- 



- 45 1 - 



m} 



10 



4 

15 4^ + 



nir 



Since they are needed for the analysis of the free energy we explicitly construct Ai & Ai 

^2 \ 3 



Ai = 
Ai = 



1 


ITIq 


12 




5 








72 



1 - 



m, 



rrin 
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C. Weak magnetic field 



1 ml 
6 



+ 



18 



6 ^10 



niQ r 



6 ^10 



r" r, 







(B.9) 
(B.IO) 
(B.ll) 



(B.12) 



(B.13) 
(B.14) 



It is possible to perform analytic computations at the limit of weak magnetic field. The definition 
of the weak magnetic field that we will adopt is that the energy scale associated to the magnetic 
field Tm is much smaller than the bare mass parameter m, namely <^ mo- In this limit to zeroth 
order in e* one has the following expansion for the profile of the fiducial embedding Xo(r) and the 
bare mass parameter mo 



Mr) = 2arcsin(^) - ^^^rt + O (r^) , 



(C.l) 



Expanding in we get 



^0 = r,--f| + 0(r^) 

9 



Mr) = ^ + + O (r^) 



(C.2) 
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which is the expression that we wiU use to calculate (3.28). Using (C.2) we can evaluate the 
expressions for all the functions of the background through (2.66), (2.71), (2.72), (2.76), (2.79) & 
(2.81) and then perform the different integrals in (3.19). After a long but straightforward calculation 
we obtain (3.28). 

Here wc give some details on obtaining the contribution from Ai, since all the rest can be calculated 
in a similar manner. Using dimcnsionless variables wc have 

JdrfV..A, ~ - ^ + + erat) log^ + OR"'). 
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